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Abstract 

The mass of the r lepton has been measured in the decay mode r — > 37r^ T using a pseudomass 
technique. The result obtained from 414 fb" 1 of data collected with the Belle detector is M T = 
(1776.61±0.13(stat.) ±0.35(sys.)) MeV/c 2 . The upper limit on the relative mass difference between 
positive and negative r leptons is \M T + — M T -\/M T < 2.8 x 1CT 4 at 90% confidence level. 

PACS numbers: 13.25.Hw, 14.40.Lb 
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Masses of quarks and leptons are fundamental parameters of the Standard Model (SM). 
High precision measurements of the mass, lifetime and the leptonic branching fractions of 
the r lepton can be used to test the lepton universality hypothesis embedded in the SM. The 
present PDG value of the r mass [l[ is dominated by the result of the BES Collaboration ^ 
and has an accuracy of about 0.3 MeV/ c 2 . The same level of accuracy in r mass measurement 
was recently reported by the KEDR Collaboration [3]. The data collected by the Belle 
experiment allow a measurement with similar accuracy to the BES and KEDR experiments 
but with different systematic uncertainties; the latter experiments analyse the cross section 
for r pair production near threshold while Belle measures the four-momenta of the visible 
r decay products at a center-of-mass (cm.) energy of yfs = 10.58 GeV. Eventually, by 
combining these high precision measurements, we will significantly improve the accuracy of 
the r mass determination. 

Separate measurements of the masses of the r + and r~ leptons in Belle, allow us to test 
the CPT theorem, which demands their equality. A similar test was previously performed 
by OPAL at LEP [4] with the result (M r + - M T -)/M T < 3.0 x 10~ 3 at 90% CL. 

To measure the r mass, we use a pseudomass technique that was first employed by the 
ARGUS collaboration [5J. This technique relies on the reconstruction of the invariant mass 
and energy of the hadronic system in hadronic r decays. The analysed variable is 



which is less than or equal to the r lepton mass. Here Mx, Ex and Px are the invariant 
mass, energy and absolute value of the momentum, respectively, of the hadronic system 
in e + e~ cm. frame, and -Ebeam is the energy of the electron (or positron) in this frame. 
In the absence of initial and final state radiation and assuming a perfect measurement of 
the four-momentum of the hadronic system, the distribution of M m ; n extends up to and 
has a sharp edge at M T . Initial (ISR) and final (FSR) state radiation as well as the finite 
momentum resolution of the detector smear this edge. We can use the edge position from a 
fit to the experimental M m [ n distribution as an estimator of the r mass, since the background 
processes in the selected t + t~ sample have a featureless M min distribution near M T . 

The analysis presented here is based on 414 fb _1 of data taken at the T(4S) resonance 
(y/s = 10.58 GeV) with the Belle detector at the KEKB asymmetric-energy e + e~ collider 
A detailed description of the Belle detector is given elsewhere (tJ • We mention here only the 
detector components essential for the present analysis. 

Charged tracks are reconstructed from hit information in a central drift chamber (CDC) 
located in a 1.5 T solenoidal magnetic field. The z axis of the detector and the solenoid 
are aligned antiparallel to the positron beam. Track trajectory coordinates of the charged 
particles near the collision point are provided by a silicon vertex detector (SVD). Photon 
detection and energy measurement are performed with a CsI(Tl) electromagnetic calorimeter 
(ECL). Identification of charged particles is based on the information from the time-of-flight 
counters (TOF) and silica aerogel Cherenkov counters (ACC). The ACC provides good 
separation between kaons and pions or muons at momenta above 1.2 GeV. The TOF system 
consists of a barrel of 128 plastic scintillation counters, and is effective in K/ir separation 
mainly for tracks with momentum below 1.2 GeV. The lower energy tracks are also identified 
using specific ionization (dE/dx) measurements in the CDC. Electrons are identified by 
combining information from the ECL, ACC, TOF and CDC [8j. The magnet return yoke 
is instrumented to form the Kl and muon detector (KLM), which detects muon tracks [9| 
and provides trigger signals. The responses from these detectors determine the likelihood 
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FIG. 1: The pseudomass distribution M m j n for the r ± — > candidates. The points with error 

bars are data and the solid line is the result of the fit with the function (2). 

Li of particle type i G {e, fi, ir, K^p}. A charged particle is identified as an electron if 
the corresponding likelihood ratio [8), P e > 0.9 or if the electron mass hypothesis has the 
highest probability. The electron efficiency for P e > 0.9 is approximately 90% for a single 
electron embedded into a hadronic event. Charged particles are identified as muons if the 
corresponding muon likelihood ratio [9| P M > 0.8. The muon detection efficiency for this 
measurement is approximately 91%. The corresponding likelihood ratio cut for kaons and 
protons is 0.8. The kaon and proton identification efficiencies are about 80%. All charged 
tracks that are not identified as an electron, muon, kaon or proton are treated as pions. The 
Kg candidates are formed from pairs of charged tracks intersecting in a secondary vertex 
more than 0.3 cm from the beam spot in the plane transverse to the beam axis; the \ 2 °f 
the vertex fit is required to be less than 11 and the it + tt~ invariant mass of the candidate is 
required to be 0.48 GeV/c 2 < M n + n - < 0.52 GeV/c 2 . The n° candidates are formed from 
pairs of photons, each with energy greater than 0.1 GeV, that satisfy the condition 0.115 
GeV/c 2 < M 77 < 0.152 GeV/c 2 . 

We select events that have one r lepton decaying leptonically into lViV T and the other 
into three charged pions and neutrino. For the entire event, we require three charged pions 
and one lepton (either muon or electron) with net charge equal to zero. The number of 
charged kaons, protons, K s mesons and 7r° mesons should be equal to zero." After applying 
all selection criteria 5.8 x 10 6 events remain. 

The M min distribution for the r — > 37ri/ T data is shown in Fig. 1. A fit was performed to 
these data with the empirical edge function 

F(x) = (P 3 + P 4 a;)arctan((x - P 1 )/P 2 ) + P 5 + P<>x, (2) 

where Pi are parameters of the fit. The fit range 1.72-1.80 GeV/c 2 is chosen. The value of 
the uncorrected r mass estimator, Pi, obtained from the fit is Pi = 1777.77 ± 0.13 MeV/c 2 . 

To obtain the value of the r mass from the r mass estimator Pi we use several Monte 
Carlo samples of r + r _ events where one r decays leptonically and the other one decays into 
three charged pions and neutrino. The KORALB generator l(j is used for the Monte Carlo 
e + e" — > t + t~ event production. 

We use three different Monte Carlo samples with input r masses equal to 1777.0 MeV/c 2 , 
1776.0 MeV/c 2 and 1776.8 MeV/c 2 for the first, second and third sample, respectively. 
The statistics of each sample is approximately equal to that of the data. The differences 
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between the fitted estimator Pi and the input r mass for these samples are Ai = (1.27 ± 
0.12) MeV/c 2 , A 2 = (1.29 ± 0.05) MeV/c 2 and A 3 = (1.06 ± 0.04) MeV/c 2 for the first, 
second and third sample, respectively. To convert the r mass estimator Pi to M T , we use 
the weighted mean and dispersion of Ai, A 2 , and A 3 to obtain the estimator correction 
A = (1.16 ± 0.14) MeV/c 2 . 

The subtraction of this value from the edge position parameter Pi in data gives M T = 
1776.61 ± 0.13(stat.) ± 0.14(MC) MeV/c 2 , where MC means the error due to limited Monte 
Carlo statistics. 

To study the systematic uncertainty due to the choice of the edge parameterization, we 
use the following alternate functions: 

Pi (:r) = (P 3 + P 4 x) x~P ^ = + ps + p (3) 

y/P 2 + (X - PI) 2 

F 2 (x) = (P 3 + P,x) — ~ l + P 5 + P 6 x (4) 

l + exp((x-Pi)/P 2 ) 

for the fit to the M min distribution. Here Pj are the parameters of the fit. 

The above procedure for r mass extraction is repeated sucessively for the data and MC 
samples with each of these functions. The extracted values for the r mass obtained with 
the functions Q and (0} are (1776.85 ± 0.13(stat.) ± 0.12(MC)) MeV/c 2 and (1776.52 ± 
0.12(stat.) ± 0.10(MC)) MeV/c 2 , respectively. 

We take for the measured value of the r mass the one obtained using function (2): 

M T = Mi = 1776.61 ± 0.13(stat.) MeV/c 2 (5) 

The square root of the variance of the obtained r masses, 0.18 MeV/c 2 is taken as systematic 
uncertainty due to the choice of the edge parameterization. This value exceeds the error 
of 0.14 MeV/c 2 , which is assigned as a systematic uncertainty due to limited Monte Carlo 
statistics. The use of a different fit range gives a much smaller shift in r mass of 0.04 MeV/c 2 , 
which we include in the systematic uncertainty. 

In this analysis, we use the beam energy calibrated using the beam-energy constrained 
mass of fully reconstructed B decays on a run-by-run basis. We estimate the uncertainty of 
the beam energy to be less than 1.5 MeV, which includes the uncertainty of B mass, tracking 
system calibration and the effect of the T(4S*) width. Using the Monte Carlo samples, we 
find that this uncertainty propagates to a systematic uncertainty on the r mass of 0.26 
MeV/c 2 under the assumption that the value 1.5 MeV is fully due to uncertainty of the 
beam energy calibration. 

As a cross-check of the result obtained from the fully reconstructed B decays we analyse 
the distribution of the variable A ME = (M(/i + /i~) — 2P bcam ) for e + e~ — > data events. 

If a systematic shift exists in the beam energy or tracking system calibration, we would 
expect some shift of the maximum of this distribution from zero. A small shift of the 
maximum of the Ame distribution from zero is due to ISR and FSR. We fit the Ame 
distribution to a sum of two Gaussians with the same central value multiplied by a cubic 
polynomial to take into account the peak asymmetry due to ISR. To check the consistency 
of this fitting procedure, we apply it to Monte Carlo e + e~ — > events with ISR and 

FSR [IH that pass through the full Belle simulation and reconstruction procedures. The 
Ame distributions for data and Monte Carlo are shown in Fig. 2 together with results of 
the fit. The reduced goodness-of-fit values x 2 /Ndf are 0.9 and 1.06 for the data and Monte 
Carlo, respectively, where Ndf = 51 is the number of degrees of freedom. 
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FIG. 2: The (M(fi + [i~) — 2Ekeam) distributions for the data (points with errors) and Monte Carlo 
(histogram without errors). The curves show the results of the fit to the data and Monte Carlo by 
the sum of two Gaussians multiplied by a cubic polynomial. 

While the resolution of the Ame variable is not well described by Monte Carlo, the peak 
position coincides for data and simulation. To estimate the systematics due to the difference 
in momentum resolution between data and Monte Carlo we included additional smearing 
of the track momenta in the Monte Carlo samples by a Gaussian with a = 1.02 • 10~ 3 p 2 
(the units of a and p are GeV/c). The consistency in Ame between data and Monte Carlo 
becomes much better after this procedure. The shift in the edge position parameter Pi is 
negligible (less than 0.02 MeV/c 2 ) and is included in total systematics. 

The difference between data and Monte Carlo peak positions obtained from the fit is 
5 Ame = 3 ± 2 MeV. This difference comes from the imperfect calibration of both the beam 
energy and tracking system. We analyse two extreme cases when the shift 5 Ame is due to 
the imperfect calibration of either 1) the beam energy or 2) the tracking system. 

For the first case, we have AE hcam = 5A M e/2 = 1-5 MeV, which is consistent with 
accuracy of the beam energy calibration obtained from the reconstruction of the exclusive B 
decays. To estimate the shift of the r mass for the second case, we construct the M m j n Monte 
Carlo distributions for an input r mass equal to 1777.0 MeV/c 2 for unmodified pion momenta 
and for momenta shifted by Ap/p = ±3/10580 = 2.8 x 10~ 4 . We obtain a mass shift in 
the range 0.10-0.15 MeV/c 2 , which is smaller than the shift observed when 5A M e includes 
the full beam energy uncertainty (0.26 MeV/c 2 ). We take this conservative assumption and 
assume a systematic uncertainty due to the combined imperfections of the beam energy and 
tracking system calibration of 0.26 MeV/c 2 . 

To estimate the systematic uncertainty due to the model dependence of the spectrum of 
the 3n system in the r decay we vary the mass and width of the ai(1260) meson in the range 
± 300 MeV/c 2 from the nominal PDG values. We find the shift in the edge position due to 
this variation to be negligible (less than 0.02 MeV/c 2 ). 

Systematic uncertainties from misidentified r decay products and from non-r + r _ events 
are negligible (less than 0.01 MeV/c 2 ), since their M min distributions show no significant 
structure in the region of the r mass. 

The list of the analysed sources of systematics is given in Table I. The final result is 
M T = (1776.61 ± 0.13(stat.) ± 0.35(sys.)) MeV/c 2 . In the analysis we assume that the 
neutrino mass is equal to zero. According to MC, a change in the neutrino mass from zero 
to 10 MeV/c 2 leads to a shift in the edge position of the pseudomass distribution by —0.1 
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TABLE I: Summary of systematic uncertainties 



Source of systematics 


a, MeV/c 2 


Beam energy and tracking system 


0.26 


Edge parameterization 


0.18 


Limited MC statistics 


0.14 


Fit range 


0.04 


Momentum resolution 


0.02 


Model of r — > 3irv T 


0.02 


Background 


0.01 


Total 


0.35 



MeV/c 2 . 

The pseudomass method allows a separate measurement of the masses of the positively 
and negatively charged r leptons. A mass difference between positive and negative r leptons 
would result in a difference in the energy between the r's produced in the e + e~ collision. 
This in principle makes the assumption E T = E heam invalid. The M min distributions for 
positive and negative r's decaying into 37tu t are shown in Fig. 3 together with the results 
of the fit. 

Good agreement is found between the distributions for r + and t~ . The mass difference 
obtained from independent fits to these distributions is M T + — M T - = (0.05 ±0.23) MeV/c 2 . 

Most sources of systematic uncertainty on r mass affect positive and negative r leptons 
equally, so that their contributions to the mass difference (and its uncertainty) cancel. One 
exception is different interactions of particles and antiparticles in the detector material. For 
example, the numbers of positive and negative triplets for the selected pions are not equal 
to each other. However, the description of this difference by the MC is reasonably accurate. 
In the data the ratio of the number of negative to number of positive triplets is 1.034 while 
in MC this ratio is equal to 1.031. To estimate a systematic shift in the mass difference 
between r + and r~ we compare the peak positions of A c — > pK~7t + and A c — > pK + 7T~ , 
D + -> 0(1O2O)tt+ and D~ -> 0(1020)71-, D+ -> 0(lO2O)vr + and Dj -> 0(1O2O)tt~ The 
average relative mass shift from the decay modes listed above is approximately 0.8 x 10~ 4 . 
This value is used as the systematic uncertainty in the relative mass difference between r + 
and t~ and corresponds to a systematic uncertainty in the mass difference of 0.14 MeV/c 2 . 

Adding the statistical and systematic errors in quadrature, we obtain M T + — M T - = 
(0.05 ± 0.27) MeV/c 2 . This result can be expressed as an upper limit on the relative mass 



difference [12j \M T + - M T -\/M T < 2.8 x 10~ 4 at 90% CL. Good agreement of the M min 
distributions for positive and negative r — > 37iu T decays shows that CPT invariance is 
respected at the present level of experimental accuracy. 

To summarize, we have measured the mass of the r lepton from the pseudomass distri- 
bution of r decays into three charged pions and a neutrino. The result is 

M T = (1776.61 ± 0.13(stat.) ± 0.35(sys.)) MeV/c 2 , 

in good agreement with the current world average and of comparable accuracy. Independent 
measurements of the positive and negative r mass are obtained to test CPT symmetry. The 
measured values are consistent and an upper limit on the relative mass difference is 

\M T+ - M T -\/M T < 2.8 x 10~ 4 at 90% CL, 
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FIG. 3: The distribution of the pseudomass M m ; n for the decay r — > 3n v T , shown separately 
for positively and negatively charged r decays. The solid points with error bars correspond to r + 
decays, while the open points with error bars are r~ decays. The solid curve is the result of the 
fit to the pseudomass distribution of r + with function (2) while the dashed one is for the r~ . 

one order of magnitude better than the previous limit from OPAL. 

We thank the KEKB group for excellent operation of the accelerator, the KEK cryo- 
genics group for efficient solenoid operations, and the KEK computer group and the Nil 
for valuable computing and Super-SINET network support. We acknowledge support from 
MEXT and JSPS (Japan); ARC and DEST (Australia); NSFC and KIP of CAS (contract 
No. 10575109 and IHEP-U-503, China); DST (India); the BK21 program of MOEHRD, 
and the CHEP SRC and BR (grant No. R01-2005-000-10089-0) programs of KOSEF (Ko- 
rea); KBN (contract No. 2P03B 01324, Poland); MIST (Russia); ARRS (Slovenia); SNSF 
(Switzerland); NSC and MOE 



[1] Review of Particle Physics, W.-M. Yao et al, J. Phys. G 33 1 (2006). 

[2] BES Collaboration, J.Z. Bai et al, Phys. Rev. D 53, 20 (1996). 

[3] KEDR Collaboration, V.V. Anashin et al, |hep-ex/0611046[ 

[4] OPAL Collaboration, G. Abbiendi et al, Phys. Lett. B 492, 23 (2000). 

[5] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 292, 221 (1992). 

[6] S.Kurokawa and E. Kikutani, Nucl. Instr. Meth. A 499, 1 (2003). 

[7] Belle Collaboration, A. Abashian et al, Nucl. Instr. Meth. A 479, 117 (2002). 

[8] Belle Collaboration, K. Hanagaki et al, Nucl. Instr. Meth. A 485, 490 (2002). 

[9] Belle Collaboration, A. Abashian et al, Nucl. Instr. Meth. A 491, 69 (2002). 

[10] S. Jadach and Z. Was, Comp. Phys. Commun. 85, 453 (1995). 

[11] S. Jadach, B.F.L. Ward, Z. Was, Comp. Phys. Commun. 130, 260 (2000). 

[12] G.J. Feldman and R.D. Cousins, Phys. Rev. D 57, 3873 (1998). 



9 



